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Abstract 



We present a method for the study of second order superhorizon perturbations in multi 
field inflationary models with non trivial kinetic terms. We utilise a change of coordinates 
in field space to separate isocurvature and adiabatic perturbations generalizing previous 
results. We also construct second order gauge invariant variables related to them. It is 
, found that with an arbitrary metric in field space the isocurvature perturbation sources the 

gravitational potential on long wavelengths even for "straight" trajectories. The potential 
decouples from the isocurvature perturbations if the background fields' trajectory is a 
' geodesic in field space. Taking nonlinear effects into account shows that, in general, the 

, two types of perturbations couple to each other. This is an outline of a possible procedure 

' to study nonlinear and non-Gaussian effects during multifield inflation. 

O ■ 1 Introduction 

^. 

\ It has been customary to say that an adiabatic, gaussian, almost scale invariant perturbation 

is a generic prediction of inflation. In the past few years though it has been realised that 
if more than one degrees of freedom are allowed to be relevant during inflation - i.e more 
dynamic scalar fields - then isocurvature perturbations can arise, possibly correlated with the 
' adiabatic ones, leading to a far richer phenomenology. For example simple statements about 

single field inflation such as the conservation of the superhorizon curvature perturbation IZ 
do not hold in multifield models pp. Therefore the need for more accurate modeling of the 
inflationary era has arisen and is actively pursued at the moment. 

So far almost all studies of inflationary perturbations have been performed using only 
linear perturbation theory. The smallness of the fluctuations in the temperature of the CMB 
certainly justify this approach. But given the accuracy of the forthcoming data it would be 
worth trying to go beyond this approximation and see if we can extract more information 
about our models by studying nonlinear effects during inflation. Nonlinearities are always 
there since gravity is a nonlinear theory. They would induce non gaussianities in the fluc- 
tuations of the cosmic microwave background which could be potentially observable by the 
PLANCK sattelite due for launch in 2007. The level of nongaussianity in standard single field 
models of inflation has been estimated in the past (e.g [21 IBIS El - ref. [01 discusses deviations 
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from gaussianity within linear theory but from a non vacuum initial state). It turns out that 
such a signal will not be detectable even by PLANCK [7j. 

Multiple scalar fields seem to have a better chance of producing an observable nongaussian 
signal |HJ the detection of which could provide evidence that more than one degrees of 
freedom were relevant during inflation. Current limits from WMAP on the deviation of 
the CMB from gaussianity come in the form of an allowed range for a widely used non- 
linearity parameter f^L related to a typs of non gaussianity. The authors of 11 find 
—58 < /nl < 134 (95% C.L.). It would be interesting to see if multiple fields can generate 
observable non gaussianity within this limit which would be observable in the future. In this 
paper an attempt is made to formulate a method for the calculation, to the lowest order, of 
the nonlinear evolution of the perturbations generated from a generic multifield inflationary 
model. This is done by extending the usual perturbation theory to second order. A future 
paper will address the issue from a different perspective jl2j . 

2 Perturbations and gauge invariance at second order 

Cosmological perturbation theory is a rather arcane subject. The reason is that in a general 
perturbed spacetime there is no privileged coordinate system with respect to which one can 
define perturbations. So perturbations can change when we change the coordinates. The 
study of general relativistic perturbations was pioneered in and studied by many authors 
since (see e.g ^1] for a comprehensive review). Let us briefiy recall a more formal presentation 
of what is usually meant when one talks of perturbations in general relativity ^SI- One 
considers a five dimensional space composed of the background spacetime Mq and, stacked 
above it, perturbed spacetimes Ai^ parametrized by the parameter e. We implicitly assume 
some sort of differentiable structure on this 5-D space such that these perturbed spacetimes 
can be considered "close" to A^o- On these spacetimes live tensor fields T. One then defines a 
vector field X, the integral curves of which are used for identifying points on with points 
on the background A^o- The choice of X is completely arbitrary and is called a choice of 
gauge. In general, any tensor T can be expanded as a taylor series 

To + SxT = <P*xeiTe) = To + e£xT |o +^e^£x£xT |o +... , (1) 
or, calling the various terms the perturbations at various orders, 

03^,(r,) =ro + 5rW + (5r(2) + ..., (2) 

where (j)'^^ is the pullback along X on the background manifold Aio of a tensor that lives 
on a perturbed spacetime Ai^ parameter distance e away from the background. Hence the 
vector field X alows us to define perturbations in a meaningful way. The choice of another 
vector field Y = {1,Y^) defines a different gauge and one finds that perturbations differ when 
defined in different gauges: 

SyT - 6xT = e£^(i)T |o +^e^i:^(i)r |o +e^£^m5^x^T + ]^e^£^(2)T |o +... , (3) 

where = e£xT |o is the linear perturbation of T in the "X gauge" and ^^^^ = Y — X, 

= [X, Y] are vector fields which lie on Mq and are independent of each other. Hence, by 
a suitable choice of ^^^^ and a gauge condition can be imposed order by order. 
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Expansion (3) also suggests a strategy for identifying gauge invariant quantities at second 
order. Observe that the first and fourth terms on th r.h.s of (3) are essentiahy the same 
(the transformations they define have the same functional form). This means that any linear 
combination f{5T^^^) of first order variables which is gauge invariant to first order will also be 
gauge invariant w.r.t that part of the second order transformation which corresponds to the 
^e'^£^(2)T |o term in (3). The remaining terms, ^e^X'^^j^^T \q +e'^£^{i)6x^T, are all composed 
of products of first order quantities. So in seeking gauge invariant combinations at second 
order we must look for appropriate quadratic terms of first order quantities that will cancel 
these quadratic terms in (3). If this can be done in a unique way then the form of a gauge 
invariant quantity at first order will dictate its form at second order. 

We will now give an explicit example of the construction of a second order gauge invariant 
variable corresponding to the well known first order gauge invariant quantity (first introduced 
in ^H]) see also [Tlj l 

7^ = V + -Y^^- (4) 

^0 

In general every quantity will be expanded in orders like in (2). For example 



^ = ^(0)+5^(l) + i5^(2) + ... (5) 

e.t.c. In particular, writing the general perturbed metric element as 

ds"^ = a2[(l + 2(l))dT'^ - 2Bidx'dT - [(1 - 2tp)-fij + 2Eij]dx'dx^] (6) 



we have 

500 



i(r)2(l + 20«+0(2) + ...) (7) 



90. = a(r)2fi?f) + i5f )+...) (8) 



9ij 



{rf [(l - 2^(1) - ^(2) + ...) + 2Eff + + ...]. (9) 



Then, from eqn. (3) one can calculate the formulae for the gauge transformations of the 
relevant quantities. For an extensive account of second order gauge transformations, explicit 
formulae and some specific examples the reader can see jl7l I18| and references therein. 

In general, formulae for perturbations at second order can be complicated and calcula- 
tions rather tedious. In this paper we will make a number of simplifying assumptions. We 
will ignore vectors (hence spatially indexed quantities are given by derivatives of scalars) and, 
mainly, we will drop terms containing more than one spatial gradients. They are expected 
to be unimportant on scales longer than the bubble radious. Although there is no rigorous 
justification for the latter approximation it is expected to capture the main affects on su- 
perhorizon scales jH]. Within such an approach, initial conditions at horizon crossing can 
be set by linear theory. Then, the long wavelength equations can be used to calculate the 
nonlinearities induced during the superhorizon evolution. The authors of ^01 adopted such 
a procedure and showed that it is possible for significant nongaussianities to be generated 
in the adiabatic mode from the long wavelength evolution in multifield inflationary models. 
Their calculation ignored metric perturbations which are included here (see next section). 
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With these approximations in mind we have the following formulae for a gauge transfor- 
mation: At first order 



^(1) 




(10) 






(11) 


V'(i) 




(12) 


E(i)ij 




(13) 






(14) 



and at second order 

0(2) = (P{2) + ^(2) + ^^(2) + 



(1) 



+ 



+ 2^o,;(2</>(,)+eo,;), 



(1) 

(15) 
(1)' 



^^(2) = V'(2) - WC(2) + ^(1) [2 (V'(i) + 27^V(i)) - (n' + 27^2) ^0^^ - neii) 



(2)' 



^(2)ii + 2^°i) [E[^)ij + 2nE^i)ij 



h{2) = 5V{2) + ^'^^2) + ill) {V"ill) + + ■ 



(16) 

(17) 
(18) 
(19) 



Note that, as mentioned before, the part of the transformations in (15) - (19) containing the 

vector field ^(2) is exactly the same as the first order case, eqn's (10) - (14). Prom the above 
we see that the variable (4) at second order transforms like 



n 



5.5^(2) +^(2) =2^5^(2) +^(2) ^ ^^_(^/ + 27^2) (^0 



+ 2(^V(i)+^(i)+27^V(i))^?i) (20) 



As expected the transformation contains only products of first order quantities. Therefore we 
seek to construct a gauge invariant quantity at second order by adding a quadratic combina- 
tion of first order quantities that will transform appropriately. By inspection we see that it 
must contain ^, (5(^' and V'' and it must not contain So we must have 



Noting that 



{All) + C^' + Db^') {Etp + GiIj' + H6ip') 



^' - ^' - n'f - Hf 



(21) 



(22) 
(23) 
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and that we must not have terms involving we see that we have 2 options. We either set 

.n 



D = C 



H = G 



which ehminates the terms involving in the transformation or take 

D = H = —n^'o 



^0 



(24) 
(25) 

(26) 
(27) 

and use the background equations of motion. In both cases we are forced to consider A = E, 
C = G and D = H and we end up with the same variable |5] 



C = G 



rri: 

n 



2 ■ 
P 

2 n'. 



n 



(2) 



^^(2) + -T^m 



+ 



.V''(i) + 2WV(i) + t5(i)^' 



n' + 2^2 _ 

which is invariant under the transformations (17) and (19). 



(28) 



3 Einstein equations for scalar fields at second order 

Consider the perturbed line element 

ds'^ = a'^[{l + 24>)dT'^ - 2Bidx'dT - [(1 - 2'il^)-fij + 2Eij]dx'dx^]. (29) 

Here it is understood that all quantities appearing are to be expanded as in (5). By inserting 
(29) into the Einstein equations with the relevant energy momentum tensor and keeping only 
linear order terms, one arrives at the well known equations of linear perturbation theory 
which can be symbolically represented as 

V6''^^g = 0. (30) 

Here I? is a set of linear differential operators and 6^^^g represents the perturbation variables. 
At second order there are two types of terms. The and terms quadratic in the 6^^^g^s. 

The later are supposed to be known from the solution of the first order problem. The form 
of the equations at second order will then be 

V5^'^g = Jmf] . (31) 

with D the same operator as in (30) and J a source term quadratic in the perturbations. 
Since the solution to the homogeneous equation T>6g = is known then we can consider the 
source terms J'HSg)'^] = J'[{S^'^^ g)'^] as known functions to second order. The solution of (31) 
will then have the form 

S^^^giv) = Jv-\y-x)Jdx, (32) 

with 'D~^{y — x) the appropriate set of Green's functions, i.e the second order perturbations 
will be determined entirely by the i7's. If 5^^^g is taken to be a gaussian random field, then 
(32) shows that 5^'^'> g is given by the square of a gaussian random field (J7 is quadratic in first 
order perturbations) and is non-gaussian. 
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3.1 Gravity sector 

In the longitudinal gauge, the linear perturbation of the Einstein tensor has the well known 
form (ignoring second order gradients) 



SlG'j = — ^({m' + n^)(t> + nct>' + 'ip'' + 2ni)')s_ 
1 



+ 



Elj + mEij 



(33) 
(34) 

(35) 



where the subscript 'L' stands for the linear part. Eij is a transverse traceless tensor (there 
are no vectors and we ignore the scalar part since it is second order in spatial derivatives). 
For the matter sector we will take a system of scalar fields (p^ with the lagrangian 



1 



which corresponds to an energy momentum tensor 



-GABdx^^d^'' - V{ip) 



(36) 



(37) 



Here, V{ip^) is an arbitrary scalar potential. The fields (p^ can be considered as coordinates 
on a field manifold with a symmetric metric Gab^- The background equations of motion for 
the scalar fields derived from (36) are 



or 



(p + 2n(p + (p (p iB(7 + aG v^F = U, 
where T^^q is the symmetric metric connection formed from Gab 

1 



^BC = 2^ {Gbd,c + GcD,B - Gbc,d) ■ 
and the linear perturbation of the energy momentum tensor is 



SlT\ = ^ 



-jGAB,cS^ 



C 



Gab ) ¥''^V'" + Gab^^'S^"" + a%c5^ 



JO 



5lT\ = ^GAB^^'d.b^'' , 



5lT' 



\Gab,c^<p'- -^Gab] ¥'^V'" + Gab'p"'^^'" - a'VcS^ 



C 



(38) 
(39) 

(40) 

(41) 
(42) 
(43) 



^The Lagrangian (36) is analogous to the lagrangian for a point particle in curved space which moves under 
the influence of a potential V. Here we have 4 parameters - the 4 spacetime coordinates - instead of one in the 
case of the point particle. 
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A perturbation Sf will be a tanjent vector on the field manifold. Given a basis in field 
space we have 

dip = Scp^eA- (44) 

If the field manifold is flat there is a prefered basis, say {ei,e2, ...} which makes the kinetic 
term in the lagrangian canonical. Even in the flat case, a basis {ga} can in general depend on 
the coordinates in field space. We will use such a basis below. Then, Gab = • ^ab- 
Here, • denotes a product in field space. If we expand ^, i/j and Stp"^ as in (5) then (5l is a 
linear operator acting at each order respectively. 

At second order wc will also have terms from the quadratic perturbation of the Einstein 
tensor and the energy momentum tensor which we denote by 620^1, and 62T^i,. Then at 
second order the Einstein equations become 

<^lG%[<^(2),V(2)] = ^dLT%[4>^2),i^(2),S^f2)]-^A (45) 

()iG0,[(/.(2),V'(2)] = ^<5LrO,[0(2), V(2),'^9'f2)]-4^' (^6) 
TThp (X 

,5iG^,[</,(2),V(2)] = ^^Lr^[<^(2),^(2),<^<pf2)]-^C^- (47) 



where 



A . (48) 



Bi - i(^2G\-^^,62T\\ (49) 



2 V S / 

^ ^ (l^^G^- - ^<52r,^ . (50) 

where, say, 52G^q is the quadratic perturbation of the Einstein tensor. To this order it will 
contain quadratic products of ijj^i-^ and which are considered known from the solution 

of the linear problem. (Please note that Bi here is different from the shift Bi of (11) - we are 
in the longitudinal gauge, and C has nothing to do with the C of the previous section). From 
(34), (42) and (46) we get that 

{ri<t> + iP') = ^Gab^^'^^"" - V-2(5,B,). (51) 



From (35), (43) and (47) we have 

^^G^^)\ = ^SlT^% - ^C (52) 

or 

{2n' + n'')(t> + net)' + < + = ^ [ {^Gab,cS^^ - 4> Gab^ 'p^'^^' 

+ Gab'p'^S^^' - a^VcS^^] - ^C. (53) 



3 EINSTEIN EQUATIONS FOR SCALAR FIELDS AT SECOND ORDER 



8 



where C = C^i. Now use the background relation 



Gab^^V^' = ^(W2-W') (54) 
Air 

and (53) becomes 

+ Gabv'"5v^' - aVc^V^^] - \c. (55) 
Note that (j) and ip appear only through the combination 

^' + n(t) = -5n. (56) 

Hence, the long wavelength sector of Einstein's equations does not contain information for 
the values of (p and "0 separately. All the terms containing the difference (p — ip are dropped 
under the long wavelength approximation. In linear theory it is known that (j) = ip m. the case 
of vanishing anisotropic stress ^^4). This need not be the case at second order. Yet, a non 
vanishing value ol (j) — ip should not matter dynamically since it does not enter the equations 
explicitly. It only appears as an initial condition for 67i provided by the short wavelength 
part of the system. But then, from equations (15) and (17) we see that we can set (p — tp = 
via a second order gauge transformation (^^^^^ = 0) by choosing 

^(2)' + 2We(2)=V'(2)-</'(2), (57) 

and keep goi = with 

C'i^2)=d42y (58) 

Of course E(^2)ij will be redefined via (18). Using (51) and the background equation of motion 
(39) we obtain in this gauge from (55) a constraint at long wavelengths for the quadratic 
parts 

{BiY + 2nBi = ^diC (59) 

Combining (45) and (47) we obtain an equation of motion for the gravitational potential at 
second order 

^p" + 6ni;' + 2(n' + 2re) = -^a^Vc5^^ + A-]-C. (60) 
V / ' 3 

3.2 Matter sector 

Equation (60) holds for an arbitrary number of scalar fields in an arbitrary parametrisation 
of the field manifold. To express the —^a^Vc^^'^' term on the r.h.s it will be useful to 
use a basis in the scalar field space which is 'adapted' to the background field trajectory, 
perturbations of which we wish to study. Such an idea was put forward in T9; in order 
to separate the entropy and adiabatic perturbations. Here we would like to give a more 
geometrical flavour which can be applied to an arbitrary number of fields and a noncanonical 
kinetic term. For simplicity we study a two field example with a diagonal metric. Some 
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results for three fields are given in appendix B. A similar approach utilising a coordinate free 
language was described in [20] 

Consider equation (51). It is a constraint equation which relates the evolution of the 
gravitational (metric) perturbations to the perturbations in the scalar fields. In general all 
the fields will be linked to the gravitational perturbations. But equation (51) holds for any 
choice of coordinates on the field manifold. Therefore if we choose a set of coordinates adapted 
to the background trajectory, i.e a set of field coordinates of which only one varies along the 
background trajectory and at the same time make Gab = ^AB on the trajectory, then the 
linear term on the r.h.s of (51) will contain only the perturbation of that field since the time 
derivatives of the rest will be zero. Call a the coordinate that varies along the background 
trajectory. It will define a basis vector 



d _ 1 d 
da a' dr 



(61) 



tangent to the background curve. Call the rest (N-1) coordinates Sj. Since we want the Sj's 
to remain constant as we move along the trajectory we demand 

(dsi, e^) = 0, i.e eg, • = 0. (62) 

We must also impose 

{dsi, eg.) = 6ij, i.e e^. • e^^ = 6ij. (63) 
The background equations of motion in these coordinates read 

a" + 2na' + cj'^r^^ + a^C^Vp = (64) 

and 

cj'^r^^ + a^G'^^Vp = 0. (65) 

The choices above mean that G^^ = 6^^ on the background trajectory in the Cq-, e^. basis. 
We also observe the following: the change of Cq- along the background trajectory - on which 
only a varies - will be 

^^Pe.e. = rie^. (66) 



where P is the covariant derivative operator in field space. Since eg- is taken to be a vector 
of unit length its variation along the trajectory will be a vector normal to it, i.e a linear 
combination of the . 's only. Therefore F^^ = 0. So the background equations of motion 
become 

a" + 2na' + a^d^^Vp = (67) 

and 

d'^r^V + a^S'^^Vp = 0. (68) 

Along the trajectory we will have 

= A^^eA (69) 



with 



We want to be a unit vector so 

1 = • = C = A^^A^ ^Gab (71) 
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from which we get that a' = \J <~p^' ^p^'G ab and hence 

A^. = , ^ (72) 

The coefficients A'^g. = dip^ /dsi which give the e^.'s in terms of the original vectors 

e,, = K^.^BA (73) 

should be chosen such that e^- • e^^ = and e^. • e^^ = 5ij, or 

GabA^^A^,, = 0, Gab^'^sA^'s, = Sij. (74) 

For > 3 there are less equations in (74) than the number of A^^^. 's needed to define the 
— 1 entropy vectors as we will now show. One has _|_ — 1 = ^^^"^^^ — 1 relations 

from (74) and one needs A'^(A^ — 1) functions to determine the e^-'s in terms of the e^'s. 
So we are left with 1 + undetermined functions. In three dimensions this is 1 which 

corresponds to our freedom to rotate the two isocurvature directions in a plane normal to e^. 

In practice it is staightforward to apply a Gram-Schmidt orthogonalization procedure to 
constuct a new orthonormal basis as follows^: We already know from eqns (69) and (72). 
Pick any other N-1 vectors from the original basis eA = which are not parallel to (as 
will be the case in general). Then one has N linearly independent vectors. By applying the 
Gram-Schmidt procedure we can construct N-1 orthonormal vectors which are also normal 
to eo- (see figure 1). These define the N-1 isocurvature (entropy) directions.^ What one ends 
up with is a relation between A"^,;- and K^a- 

A^s. = {a\) . (75) 

The only thing that still needs to be specified in order for this new coordinate system to 
be defined in a region around the trajectory are derivatives w.r.t the Sj's. They will be 
determined by demanding that the isocurvature part of the basis, e^. = d/dsi, commutes 
with 

[es,,e^]=0. (76) 
By acting the above commutator on the functions (p^ we get 

dsi a'dr a'dA^dr ^^'> 

with defined in (75). The derivatives dsiA^g. are of course defined via (77) and (75). 
Equation (77) can be used to compute any derivative in the isocurvature directions. 
Let us now focus for simplicity on a system with two fields and a diagonal metric 



Gab={ " : . (78) 




■^Many thanks to Christopher Gordon for a stimulating discussion on this 

■^For an apphcation of a Gram-Schmidt orthogonahsation to assisted inflation models see |21| 
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(p3 




s = constant 

background trajectory 



Figure 1: The new basis in N-dimensional field space. Here, a is a coordinate that varies 
along the background trajectory. It defines a basis vector e^- tangent to the curve. There are 
— 1 coordinates Si which are constant along the curve and define — 1 coordinate directions 
normal to it. The basis vectors along these directions e^. are taken to be orthonormal and 
define the A*" — 1 isocurvature perturbations. 



We take s = on the trajectory and we write (72) as 



A^, = U{s) , '^^ (79) 



y'a((^i')2 + /?(v.2')^ 



t^\ = W{s) , (80) 

n2'^2 



y'a((^i')2 + /3((^2')s 



where we have parametrised the dependence on s by the functions C/(s) and W{^s). On the 
trajectory f7(0) = VF(0) = 1. We take the unnormalised vector normal to the velocity to be 

6s = 62 - (e, • e2)e, = (l - /?(A^a)^) 63 - /^A^A^^ei. (81) 
Therefore, the isocurvature basis vector is 



-A^ei + ./-AVe2, (82) 
from which we read off the transformation coefficients in the s direction 

K\ = = (83) 

V a VP 

The metric in these new coordinates reads 

C = Gss = a (k\) ' + (3 (k\) ' , (84) 

Gsa = 0, (85) 

which reduces to unit diagonal on the background trajectory, s = 0, as expected. 

In order to specify the new coordinate system completely, we need to know the functions 
U and W appearing in (79) and (80), or, equivalently, to be able to calculate a derivative 
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of arbitrary order with respect to s on the background trajectory, i.e at s = 0. This can be 
achieved via (77) which in this case reads 

and 

With these, any derivative in the a, s coordinates can be calculated in terms of functions of 
r, known along the background trajectory. 

We can now derive equations for the gravitational potential and the isocurvature pertur- 
bation at second order. The first term on the r.h.s of (60) is written in the a, s coordinates 



ml mt ' m, 



p 



Prom (67) 

and from (51) 



a\^Sa = + (r'Sa, (89) 



2 2 

Tnz, , „ rat 



a' 5a = -P{n^ + + (diBi) . (90) 

47r 47r 

So from (60) we can get an equation for the gravitational potential to second order 

+ A- F' -2(^^ +3Tij F, (91) 

with 

F = V-^diBi. (92) 
and we have used that on long wavelengths we have from (59) 

F' + 2nF = \c. (93) 

We also need an equation for the entropy perturbation 5s. We can get it by perturbing (38) 
expressed in the cr, s coordinates: 

5s" + 2H5s' - 2i,{a'fTl„ + 2a' [sa'V^^ + 5s'V's^ 

+ {a'f [ssT^^^s + 6aT%^,^^ + a^G'^ ^^5v^Vp + a''V^^^5^^ = r. (94) 

Expressing 5a and 5a' in terms of ip and ip' we finally get 

5s" + 2n5s' + \a^Vss - 6 [a'f + [a'f ^sTl^ = r. (95) 
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The term appearing on the r.h.s of (95) is a second order source and it is calculated in 
the appendix. Equations (91) and (95) are generalisations to second order and a diagonal 
but otherwise arbitrary field metric of the equations for the evolution of long wavelength 
curvature and isocurvature perturbations given in ^Hl- The above procedure could of course 
be applied to an arbitrary number of fields and the most general metric. 

A nice result proved in _19' is that when the background trajectory is straight in field 
space then isocurvature perturbations do not source the adiabatic one or, equivalently, the 
gravitational potential perturbation. This is true only if the metric is flat, Gab = ^AB- 
If the field space has a nontrivial metric then the equivalent statement would be that this 
decoupling occurs if the background fields follow a geodesic in field space. By geodesies we 
mean the curves that the background fields would follow if there was no potential term in 
the lagrangian (36). We will now prove this for the two field case but the argument can be 
generalised in a straightforward manner. The geodesic equation in field space would read 

+ 2n^^' + ^'"^'"T^BC = 0. (96) 

In the o", s coordinates these equations read 

a" + 2n(j' = 

{a'fr,„ = 0. (97) 

Hence, such trajectories correspond to F'^o-o- = and from (91) we see that, indeed, the isocur- 
vature perturbations do not source the metric perturbation ^. In contrast, for a "straight" 
trajectory T^o-o- 7^ 0. Indeed 

— 2 (^^S(T,(T Gcra,s) 

= -\Gaa,s (98) 

where we have used that, by construction, = (^^^ along the trajectory. Now G^a = 

A^^A^^Gab so 



Gaa.s = '^A'^^A^ sGab + A'^cr A^cr A*" sGab, 



c 



= 2GAB^^a-,^A^, + A^„A^,A''sGAB,c (99) 
a' OT 

from eq (77). In general the r.h.s of (99) will be different from zero. Hence T^o-o- 7^ and 
isocurvature perturbations can source the adiabatic one. 

We conclude that nontrivial kinetic terms can lead to nontrivial couplings between the 
adiabatic and the isocurvature perturbations which may not be supressed on superhorizon 
scales, in contrast to what happens in models with flat field metrics ^ . Such scalar Lagrangians 
appear in the effective actions of various candidate fundamental theories. Therefore it would 
be worth studying wether the interplay of isocurvature and adiabatic perturbations in such 
models has any interesting phenomenological consequences. Of course, non linear evolution 
will couple the two types of perturbations anyway. 

*r^CTcT measures the rate of change of the vector along the trajectory in the direction of s. Hence it is 
true that a trajectory that bends generates a coupling between isocurvature and adiabatic perturbations as 
suggested in 19 • 

particular example for a two field model with a non flat metric where such a conclusion is reached was 
studied in |22 
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3.3 Superhorizon Gauge Invariant Variables 

Having found solutions to the perturbation equations at second order in a particular gauge, 
we can always construct gauge invariant quantities in terms of the new field coordinates a 
and Si. At first order the gauge invariant curvature perurbation TZ is can be written as 



(1) 



1~C 



a' 



(100) 



Since s' = 0, 5si is gauge invariant at first order. At second order, according to eqn. (19) 

5s = (5s + 2e['i)<5s/(i) (101) 
so a corresponding gauge invariant quantity is easily seen to be 



6 Si 



2— — (^Si(i). 



CJ' 



(102) 



Similarly we can construct a gauge invariant curvature perturbation at second order similar 
to (28) 



7^ 



(2) 



V'(2) + -,^^{2) 



+ 



V''n)+2W^(i) + |^5(i)a' 



W + 2^2 _ 



(103) 



The variables (102) and (103) are second order gauge invariant variables (at least on super- 
horizon scales) which can be used to study isocurvature and adiabatic perurbations in the 
mild nonlinear regime. 



4 Summary 

We have touched upon the issue of studying second order perturbations in multifield inflation- 
ary models and defined gauge invariant variables - equations (102) and (103) - on supehorizon 
scales also to second order. The latter can be constructed given the solution to lowest non 
linear order in a given gauge. We have presented a more geometrical method for the splitting 
of isocurvature and adiabatic perturbations and applied it to a two field model with a diago- 
nal but otherwise arbitrary metric. We showed that in this case naive "straight" trajectories 
do not lead to the decoupling of adiabatic from isocurvature perturbations. We identified 
the type of curves for which this happens. A perturbative approach by which one can study 
the nonlinear evolution of perturbations in such models to second order was suggested. The 
resulting equations have the same form as the first order linear ones but with new terms 
appearing on the right hand side. These new terms are quadratic in the first order pertur- 
bations and therefore they can be considered as known "sources" from the solution of the 
first order problem. Such a formalism can be used to calculate the amount of nongaussianity 
produced in such models by treating the perturbations as gaussian stochastic fields when they 
become superhorizon and then study their non linear evolution from that point. The resulting 
non gaussianity will be of the \^ type since we are considering quadratic products of gaussian 
fields. We have implicitly assumed a smoothing on scales larger than the horizon and dropped 
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second order spatial gradients. Although straightforward to calculate, the resulting source 
terms, in principle known, are quite complicated. We will return to the issue of calculating 
non linear evolution in multifield inflationary models with a different approach in a future 
publication T2 . 

Acknowledgements: Many thanks to Carsten van de Bruck and Christopher Gordon 
for useful discussions and comments, Marco Bruni for bringing to my attention other works 
related to second order perturbation theory, and especially Paul Shellard for discussions and 
generous support. 



A Second order sources 

In this appendix we give formulae nessecary for the calculation of the terms A, C, F, J* 
defined in the text. These are the terms containing products of first order perturbations. 
We write the metric as 

where /iQi = 0. Indices will be raised and lowered with the background metric. Then the 
perturbation of the contravariant metric tensor will be 



and in particular 



The perturbation of the metric determinant is 



9' 



5°° = \- -Khoo + ^/loo/ioo (106) 
g^' = (107) 
3^^ = -% - -ihij - \hikhkj. (108) 



^=a^[l + ]^h+ - ^h^ph^^j . (109) 

The perturbation of the Riemann tensor to first and second order can be found to be f j23j 
p:965) 

-^Mi'^ = 2 {^^\^^'' ~ Vi'Iq" + ^omIi^" + (110) 

and 



2 

+ 



fill) 



where a vertical bar denotes a covariant derivative w.r.t. the background. Then 

RW^'^ = - (112) 

and 

i?(2)/^^ = g(^)f^-R(^} - /i'^^i?^)^, + h^xh^'^Ri'll (113) 
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Using the following perturbed metric tensor (longitudinal gauge) 

hoQ = 2a^il} 
hij = 2a'^'^5ij — 2c?Eij 



we get, after a rather tedious calculation 
2 



6WVV -Q{n' + 2n^) i;^ - Eki (UE'^i + 2H''Eki) + ^E'^iE'^i 



di^ (^' - m^) + 2Eki {diE'^i - diEi,) + E',idiEki + 2Eudii;' - 2di'il^E'^ 



(114) 

(115) 
(116) 

(117) 
, (118) 



+ ^Eij + 2nip' + 2{m' + m^)^) 

~ o'^^'ik^kj - Ti^E'^iEki^ii - -{H' + 2'H^)EikEkj 



The Einstein tensor will be 



so we find 



^(2)0 . ^ ^(2)0 . 
Q{2)i. ^ j^(2)i. _ 1 . (^^(2)0^ _ ^(2)i . j ^ 



(119) 

(120) 
(121) 
(122) 



(123) 



+ E'kAEki + 2EiM - 2diijEii 



(124) 



+ "^'^^'ki^kj - 'HE'i^iEkidij + {H' + 2Tt^)EkiEkj 

- 'jp^^ ~ '^T-C')EkiEkidij - -E'kiE'^.iSij 

- 2Eij (V" + 2n^^J' + 2(4W' + dU' 



(125) 
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For the energy momentum tensor 



the second order perturbation for two fields is given by 



{Saf - 2a' (KJ SsSa - {a')' dsKJs 



/\2 



+ 2iP^ {a'f + 2 {a'^T'^^ipSs - 2T'^y6sda' - 2a'ipda' 



+ VssSs'^ + 2Vsa5sSa + V^^Sa'^ 



t(2)o . = [5s'di5s + {5a' - 2a'T'^Js - 2iPa') diSa] , 



(126) 



(127) 

(128) 



{Sa'f - 2a' {FU ^sda - {a'Y dsKJs 



/^2 



+ 2i;^ {a')^ + 2 {a'Y F'^^il^Ss - 2T'^yds5a' - 2a'ip5a' 



/^2 



+ 5) 



(129) 



where wc have obviously dropped second order spatial derivatives. We also need to compute 
the quadratic part of the wave equation for the system of scalar fields 



Up to second order spatial gradients we get 

A 

+ «V'(nj' 



js = -4^(a'^6s')+4a^a'fT^^^{i;' + K,^ds) 



?,^5a + 2^-^ + 2^-^+AKJs5a 



a' 



2 



(130) 



(131) 



-26s'^ - 3 {a'fi>5s + 2a'5s6a' + 2a'6sSs' 

+ a'{a'y6a 

- 2aV {VssSs + VsJa) + a" {VsssSs^ + 2VssaSsda + VsaaSa^ 
With these formulae we can compute all the second order terms defined in the text. 

B Isocurvature-adiabatic split with three fields 

In this appendix we present the isocurvature-adiabatic split for three fields (p^, ip^, ^p^ with 
an arbitrary metric Gab- This should illustrate the general case. Now we will have two 
isocurvature directions. Take the set {e^, 62, 63} as the starting point for constructing the 
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orthonormal vectors. Form e^j which hes in the plane spanned by and 62- First form the 
unnormaUsed 

= e2 - (62 • ea)ea- 



One finds 
with 



-S2 



S2 



A' 
A' 



S2 



where 
Then 

So, we have 



1 - A2aA^a = Ai^A^a + As^A^,, 
-A2aA^a, 



Aao- = Gab A a- 

A^ 



-'52 



VA^^Ab^ 
A^ 



S2 



Repeating the procedure we form 

= 63 - (63 • eCT)e(^ - (§3 • e^Jesj 



to get 



A^3 = 



S3 



A3 



S3 



-Aso-A^o- - AasjA^sj 

— A3(tA^0- — A352A^52 
1 - Asn-A^^ - A3.<=oA3 



So we have 



A" 



A^^„ 



As before, we demand for the e^.'s 



i^si, ea] = 



BS3 



from which we get by acting on ip^ 



— — A^.„ = 



a' dr 



a'dT ds3 



d 

OS2 

a- 



(132) 
(133) 



(134) 
(135) 
(136) 

(137) 
(138) 



(139) 

(140) 

(141) 

(142) 
(143) 



Knowing the A coefficients we can calculate the metric in the new basis, where the index 
A runs over {a, si, 82}- Using equations (142) and (143) we can then calculate any partial 
derivative of A^cr in terms of time derivatives along the background trajectory. Hence we can 
calculate the connection and perturbations of arbitrary order in the new {e^, e^^, 652} basis. 
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